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We consider a bosonic Josephson junction made of A*" ultracold and dilute atoms confined by a 
quasi one-dimensional double-well potential within the two-site Bose-Hubbard model framework. 
The behaviour of the system is investigated at zero temperature by varying the inter-atomic in- 
teraction from the strongly attractive regime to the repulsive one. We show that the ground-state 
exhibits a crossover from a macroscopic Schrodinger-cat state to a separable Fock state through 
an atomic coherent regime. By diagonalizing the Bose-Hubbard Hamiltonian we characterize the 
emergence of the mascroscopic cat states by calculating the Fisher information F, the coherence by 
means of the visibility a of the interference fringes in the momentum distribution, and the quantum 
correlations by using the entanglement entropy S. Both Fisher information and visibility are shown 
to be related to the ground state energy by employing the Hellmann-Feynman theorem. This re- 
sult, together with a perturbative calculation of the ground-state energy, makes possible to obtain 
simple analytical formulas for F and a over a range of interactions, in excellent agreement with the 
exact diagonalization of the Bose-Hubbard Hamiltonian. In the attractive regime the entanglement 
entropy attains values very close to its upper limit for a specific interaction strength lying in the 
region where coherence is lost and self trapping sets in. 

PACS numbers: 03.75.Ss,03.75.Hh,64.75.-l-g 



I. INTRODUCTION 



Ultracold and dilute alkali-metal atoms confined in a 
quasi one-dimensional (ID) double-well potential_[l| arc 
the ideal system to study the Josephson effect and, 
more generally, the formation of macroscopic coherent 
states [^-01 and macroscopic Schrodinger-cat states 
[T^ . The observed coherent dynamics of the atomic cloud 
in the double- well potential (bosonic Josephson junction) 
1, [l^ is efficiently described by Josephson equations 
3 1, and their extensions [l3 - [l6j . The Josephson equa- 
tions are valid within the weak inter-atomic interaction 
regime, where scmiclassical approximations can be per- 
formed and the system considered is in a coherent state 
[l7| . When the bosons are repulsively interacting, the 
crossover from a coherent state (superfluid-like regime) 
to a separable Fock state (Mott-like regime) takes place 
by increasing the coupling strength between the atoms 
[^-Q . This crossover could be detected by measuring the 
visibility a of the interference fringes in the momentum 
distribution of the bosonic cloud ]3, 0] ■ ^^e case of 
attractive coupling, the Josephson equations predict the 
spontaneous symmetry-breaking above a critical strength 
S EBl J while the two-site Bose-Hubbard model predicts 
the formation of a macroscopic Schrodinger-cat state @- 
[l^ . Finally, when the attraction between the bosons is 
sufficiently strong the collapse of the cloud should occur 




The goal of the present paper is to investigate the co- 
herence and entanglement of the N boson ground-state, 
with particular attention to the formation of macroscopic 
Schrodingcr-cat states, by using as theoretical tool the 
two-site Bose-Hubbard (BH) Hamiltonian. Both the ex- 
perimental implementation of two-mode quantum entan- 
gled states and the two-site BH model are the ideal in- 
struments to address the problem of the precision of the 
difference number and phase measurements p^ . [20j . In 
our work, we diagonalize the two-site BH Hamiltonian 
and evaluate the coherence visibility a, the Fisher infor- 
mation F, and the entanglement entropy S of the ground 
state. We calculate the Fisher information and the co- 
herence visibility by employing the Hellmann-Feynman 
theorem (HFT) (see, for example, Ref. [2lj) as well. By 
jointly using this theorem and a perturbative approach 
for calculating the ground-state energy, we obtain simple 
analytical formulas for F and a over range of interactions 
where the perturbative calculation is applicable. This 
represents one of the main novelties of the present paper. 
It is quite interesting to compare the results obtained 
by applying the HFT and the ones deriving from the di- 
agonalization of the BH Hamiltonian. This comparison 
shows a good agreement between the predictions of the 
two approaches, in particular for the Fisher information. 

The aforementioned quantities, i.e. the Fisher infor- 
mation, the coherence visibility and the entanglement en- 
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tropy are analyzed by widely exploring the atom-atom in- 
teraction range, i.e. from strongly attractive to strongly 
repulsive interactions. This study is - in particular, the 
analysis of the entanglement entropy - one of the crucial 
points of our paper. From the above analysis we find that 
the presence of a macroscopic Schrodinger-cat state [22| 
in the double-well, which can be detected by repeated 
measurements of the distribution of particles in the two 
wells d,!!!, corresponds to the maximum of the Fisher in- 
formation F. The maximum of the entanglement entropy 
S is obtained for attractive interactions close to the on- 
set of coherence loss, which corresponds to the crossover 
to the self-trapped state. This maximum value is very 
close to its maximum achievable value, i.e. log2(A^ + 1), 
and greater than the one of a macroscopic cat state. For 
sufficiently strong attractions the coherence visibility a 
is quite small but it increases as the interaction strength 
becomes sufficiently weak. In addition, we find that, by 
strongly increasing the repulsive interaction, the coher- 
ence visibility, the Fisher information F and the entan- 
glement entropy S go to zero. In this strongly-repulsive 
regime the ground state becomes a separable Fock state 
with N/2 particles in each potential well. 

To make complete our analysis, we exploit the known 
fact that the Josepshson equations for a bosonic junction 
can be derived from the two-site Bose-Hubbard Hamilto- 
nian by using the quasi-classical coherent state; this state 
is characterized by non- vanishing quantum number fiuc- 
tuations. We compare the A^-boson ground-state energy 
obtained by diagonalizing the BH Hamiltonian with the 
expectation value of the BH Hamiltonian evaluated with 
respect to the quasi-classical coherent state. In this way, 
we show the link between the quantum description and 
the quasi-classical one of the bosonic Josephson junction. 



II. THE MODEL HAMILTONIAN 

We consider a dilute and ultracold atomic gas made 
of N identical bosons of mass m confined by a trapping 
potential Vtrapi'^c)- This potential can be realized by the 
superposition of an isotropic harmonic confinement in the 
the transverse radial plane and a double-well potential 
Vdw{x) in the axial direction x. Then, Vtrap(r) is given 
by Vtrap{r) = Vdw{x) + muj]_{y'^ + z^)/2, where wj. is 
the trapping frequency in the radial plane. 

We suppose that the system is quasi one-dimensional 
(ID) due to a strong transverse radial harmonic confine- 
ment. In particular, the transverse energy Tiloa^ is much 
larger than the characteristic energy of bosons in the ax- 
ial direction. By assuming that the two wells are symmet- 
ric, the microscopic dynamics of our system is described 
by the effective two-sites Bose-Hubbard (BH) Hamilto- 
nian uM 



where afc,aj, (fc = L,R) are bosonic operators satisfying 



H 



J {a\aii + aJjOi) + ^ (nL(nL - 1) + nR{nit - 1)) 



U 



the algebra [hk , a\] 



Oki; rik 



alak is the number of 



particles in the fcth well; U is the boson-boson interaction 
amplitude; J is the tunnel matrix element between the 
two wells. Note that the sign of the hopping amplitude 
can be changed by gauge transformation: in the following 
we will refer to J as a positive quantity. We observe that 
the total number operator 



N ^ fiL + hji 



(2) 



commutes with the Hamiltonian ([T]), so that the total 
number of particles is conserved. The macroscopic pa- 
rameters in the Hamiltonian ([T]) are explicitly related to 
the atom-atom coupling constant g = ATrh^as/m - as 
being the s-wave scattering length - and to the other mi- 
croscopic parameters, i.e. the atomic mass m and the 
frequency ujj_ of the harmonic trap (see Ref. [l^ ) . The 
on-site interaction amplitude U is positive (negative) if 
tts is positive (negative), so that it may be changed by 
Feshbach resonance. The hopping amplitude J is equal to 
(ei — eo) /2, where eo and ei are the ground-state and first 
excited state energies of a single boson in the double- well 

ii. 



HI. COHERENT, FOCK AND CAT STATES 

At fixed number N of bosons, the ground-state of the 
system depends on the adimensional parameter ^ = U/ J. 
Let us recall what is known in the literature. 

• In the case of zero interaction {U = 0) the ground- 
state of the BH Hamiltonian with a fixed total num- 
ber N of bosons is the atomic coherent state (ACS) 



\ACS) 



1 



N 



|0,0) 



(3) 



where |0, 0) = \Q) r, that is the tensor product 
between the vacuum of the o per ator ol - |0)l - and 
the vacuum of and an, |0)ij |24l. [25j . 

• In the case of repulsive interaction {U > 0) the 
ground-state of the BH Hamiltonian can be found 
into two very different regimes: the superfluid-like 
regime and the Mott-like regime (343 • The transi- 
tion from one regime to the other depends on the 
parameter ( = U/ J > 0. li ( <SC. N, the system 
is in a quasiclassical superfiuid regime character- 
ized by the suppression of phase fluctuations: the 
ground-state is close to \ACS). If ( ^ N, the 
BH Hamiltonian ([T]) yields the Mott-insulator like 
regime where number fluctuations arc frozen: the 
ground-state is incoherent [l-Q, and close to the 
separable Fock state 



(1) 



N N 

\FOCK)^\-,-) 



(4) 
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• In the case of attractive interaction {U < 0) the 
ground-state of the BH Hamiltonian can be close to 
the atomic coherent state \ACS) or in a entangled 
superposition of macroscopic states, i.e. macro- 
scopic Schrodinger cat states [Ml. By changing 
the parameter ( = U/ J < the ground-state of the 
system evolves towards the following macroscopic 
superposition state 



\CAT) 



1 



V2 



(|^,0)-f |0,iV)) 



(5) 



which is the linear combination of the states with 
all particles in the left or in the right well. This 
entangled macroscopic superposition state is also 
known as ""NOON state" or "macroscopic cat 
state" . The lower is C the more likely the system 
is expected be in the macroscopic Schrodinger-cat 
state [Ml- 

A very interesting issue to address is the study of the 
coherence and the entanglement in a bosonic Josephson 
junction by varying the ratio U/J from negative to pos- 
itive values keeping fixed the number N of atoms. To 
accomplish this task we solve numerically the following 
eigenproblem 



H\E,) = E, 



(6) 



for a fixed number TV of bosons. In this case the Hamil- 
tonian H can be represented by a {N + 1) x [N + 1) 
matrix in the Fock basis \i^N — i) = \i)L®\N — i) b. with 
i = 0, N . For each eigenvalue Ej, with j = 0, 1, N , 
the corresponding eigenstate \Ej) will be of the form 



\E, 



N 

E 

4=0 



i,N - 



(7) 



In Fig. [T]we plot the coefficients |c,|°''p of the repulsive 
[U > 0) ground-state \Eo) for different values of the pa- 
rameter ( ^ U/ J. In correspondence to each value of 
we have studied the ground-state for N = 30 (solid line) 
and N — 300 (dashed line). Let us focus on the solid line. 

The figure shows that the width of the peak of {cf^^ p be- 
comes narrower by increasing C. The hardening of the 
localization corresponds to a reduction of the coherence. 
This is related to the fact that the ground-state \Eo) is 
becoming the separable Fock state \FOCK). The same 
kind of behavior is observed if the number of particles 
N is increased. From the plots reported in Fig. [U we 
can see that, given the greater is N the stronger is the 
localization. 

In Fig. [2] we plot the coefficients |Cj-''''p of the attrac- 
tive {U < 0) ground-state \Eo) for different values of the 
parameter ( ~ U/J. In correspondence to each value of 
C, we have analyzed the ground-state for A'' = 20 (solid 
line) and = 30 (dashed line). Once fixed the num- 
ber of bosons A^, we see that by increasing ( a crossover 
occurs. It exists a value of say Ccr (Ccr — —0.116 for 
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FIG. 1: (Color online). Coefficients |c''''p of the repulsive 
(U > 0) ground-state \Eo) (see Eq. (O) of the bosonic junc- 
tion as a function of i/N. \ci,of are plotted for four values 
of the adimensional parameter ( = U/J > 0. Note that 
X/i^o P = 1 ^^'^ the total number of coefficients is N + 1. 
The coefficients |cf' ^ and i/N are adimensional quantities. 
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FIG. 2: (Color online). Coefficients Ic'^'p of the attractive 
(U < 0) ground-state \Eo) (see Eq. ((TJ) of the bosonic junc- 
tion as a function of i/N. \ci,of are plotted for four values 
of the adimensional parameter = U / J < 0. The coefficients 
|c^°'P and i/N are adimensional quantities. 



A^ 20 and Ccr ^ -0.076 for A^ 30), for which a valley 
appears in the middle of {cf' p. This is a signature of the 
fact that the system is losing coherence and the junction 
is approaching the self-trapping regime. The above men- 
tioned crossover takes place from a ground-state which 
has the maximal probability at i ~ N/2 (single peak 
centered around i ~ 0) for C = 0, to a ground-state 
which has the same maximal probability at i = n and 
i = N — n (two separated peaks symmetric with respect 
to i = 0) with n approaching as C decreases. Thus, 
for a sufficiently small C (i-e. C < ^ind |C| large) the 
ground-state \Eq) has the maximal probability to be in 
the state in which the bosonic population in each well 
fluctuates around A^/2. For a sufficiently large instead. 
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the ground state \Eo) has the maximal probabihty to be 
in a "macroscopic cat state" : a Hnear superposition of 
states with macroscopically different occupations. From 
the plots shown in Fig. [21 wc can see that the larger is 
A^, the larger is the attraction for which the crossover to 
the "macroscopic cat state" occurs. 

Finally, we can summarize that Fig. [T]and Fig. [5]show 
the crossover from the "macroscopic cat state" \CAT) for 
C < to the "separable Fock state" \FOCK) for C > 0, 
crossing the atomic coherent state \ACS) at C = 0. 

It is quite important to note that the ground-state \Eo) 
has definite parity, namely the distribution of the coeffi- 
cients cf ^ is even with respect to i/N = 1/2. As conse- 
quence, we have that {ul) = N/2 ~ {fiR). 

To conclude this section we want to stress the role 
played by the gap, Aii^, between the ground-state and the 
first excited state of the Hamiltonian ([1]). In fact, when 
this gap approaches zero the even and the odd states be- 
come quasi-degenerate allowing parity symmetry break- 
ing. Let us focus on the attractive interactions, U < i). 
In the classical limit, U — >■ — oo, the ground state is dou- 
bly degenerate comprising the two Fock states |iV, 0) and 
|0, N). The degeneracy is lifted by quantum fluctuations 
and a gap between the even and the odd linear combi- 
nation of the two Fock states sets in. By standard semi- 
classical analysis, see Ref. H^, we can estimate the gap 
Ai? between the two quasi-degenerate ground states in 
the \U\N oo limit: 



\U\^N{N + 2) 



where ^ is a positive quantity given by 



i = iog( 



\U\N , 
J ' 



(8) 



(9) 



This asymptotic expression holds for sufficiently large 
values of the product UN, as shown in Fig. [31 where 
IS.E is plotted as a function of N for U/ J = —0.15. We 
can see that the diagonalization of the Hamiltonian ([T]) 
predicts a gap closing in according a law - dashed line 
- which is well represented by the one given by Eq. ([5]), 
solid line of Fig. [31 



IV. COHERENCE AND ENTANGLEMENT AT 
ZERO TEMPERATURE 

In this section we shall quantitatively characterize the 
coherence and the entanglement properties of the ground- 
state, by devoting particular attention to the emergence 
of the "macroscopic cat state" . Let us start with this last 
issue. To this end, we recall the definition of the quantum 
version of the Fisher information Fqpi [22l. [27l. [28| 



Fqfi = {AhL,Rf 



where 



(10) 



(11) 




FIG. 3: (Color online). On the vertical axis the energetic 
gap Ai5 at U / J = —0.15. On the horizontal axis the number 
of bosons N . The solid line is obtained with Eq. The 
circles on the dashed line are the points obtained from the 
diagonalization of the Hamiltonian Note that the vertical 
axis is in logarithmic scale. Ai5 is in units of J. 



It is worth observing that the above definition of the 
Fisher information, that we use hereafter, only holds for 
pure states (i.e. only for not mixed states). The expec- 
tation values at the right hand side of Eqs. (|TU| and 
Eqs. (|lip are taken with respect to the ground state 
\Eq) which, having definite parity, satisfies the condition 



{til) = {nR) 
given by: 



In terms of the coefficients c, 



(0) 



F, 



QFI 



N 



QFI 



1=0 



(12) 



It is convenient to normalize Fqpj at its maximum value 
N'^ by defining the Fisher information F as 



F = 



F, 



QFI 

7V2 



(13) 



F gives the width of the distribution |c^''''|2 centered 
at i/N — 1/2. As expected, F =1 holds for the "macro- 
scopic cat state" \CAT), and _F = for the "separable 
Fock state" \FOCK). 

We have studied F as a function of the parameter 
C = U / J . The results are shown in Fig. [H for iV = 20 
(solid line) and TV = 30 (dashed line). The Fisher Infor- 
mation is a monotonic decreasing function of Q both 
for attractive and repulsive interactions. On the attrac- 
tive side - left panel of Fig. [H- it exists a value of say 
Cf (Cf ^ -0.135 for TV = 20 and ~ -0.08 for iV = 30) 
for which the second derivative of F is zero. 

It is quite interesting, now, to observe that by using 
the Hellmann-Feynman theorem (HFT) - see, for exam- 
ple, Ref. [U - a relation can be established between the 
Fisher information F and the partial derivative of the 

dE 

ground-state energy Eq with respect to C/, "q^j- ^c- 
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FIG. 4: (Color online). Fisher information F of the ground- 
state l-Eo) of the bosonic junction as a function of the param- 
eter = U/J. Left panel: attractive bosons {U < 0). Right 
panel: repulsive bosons {U > 0). Solid line: N = 20. Dashed 
line: A'^ = 30. Note that the horizontal axis of the right panel 
is in logarithmic scale. The Fisher information F and ( are 
adimensional quantities. 



cording to the HFT, we have that 
dEo _ dH 



(14) 



The aforementioned relation is given by the folloviring for- 
mula: 



F : 



4 dEo 



2 



(15) 



In the limit U/J — > — cx3, the hopping operator in Hamil- 
tonian ^ can be treated within perturbation theory. To 

the zero order, J = 0, the energy e'^q^ of the ground- 

UN 

state of the on-site interaction operator is —^{N — 1). 

This level is twofold degenerate, and the two states cor- 
responding to this energy satisfying (rii) = N/2 = {fiji), 

i.e. i=(|Ar,0) + |0,iV)) and 1={\N,0) - \0, N)), have 

opposite parity. 

At the first non vanishing order the ground state en- 
ergy is given by [29| : 



J2 N 



U N -1 



(16) 



By substituting this energy into the expression ()15p . we 
obtain 



F = 1 



N{N-1)C^ 



(17) 



which approaches to 1 as C — > — oo. 

Also in the limit U/J^ +oo, we can treat the hop- 
ping operator of the Hamiltonian ([T]) as a perturbation 



to the J = BH Hamiltonian ([T]). To the lowest order, 

the ground state energy Eq^^ of the on-site interaction 
UN N 

operator is equal to ^ (— — 1), which corresponds to 

the separable Fock state (j^. At the first non vanishing 
order we get 



(0) J^N ,N . 

Eo^eT-^{^ + i) 



(18) 



By substituting this energy into the expression (jlSp . we 
obtain 



F 



2{N + 2) 
NC 



(19) 



which goes to zero as — >■ -foo. 

In the opposite limit, [/ — )■ 0, the on-site interaction 
term can be considered as a perturbation in the Hamilto- 
nian ([1]). For U ~ 0, the ground-state energy Eq^^ of the 
hopping operator is equal to —N J which corresponds to 
the atomic coherent state (jS]). To the first order in per- 
turbation theory, we obtain 



En ~ E, 



(0) , UN{N-l) 



(20) 



By substituting this energy into the expression (jlSp . we 
obtain 



-4 



(21) 



which tends to when N tends to infinity. In Fig. ([Sj we 
show the comparison between the Fisher information F 
obtained by diagonalizing the Hamiltonian ([1]), dashed 
line, and that one obtained by using Eq. (jl7p . attrac- 
tive bosons, and Eq. ([T9| . repulsive bosons, see the dot- 
dashed line. 

In cold atom physics it is customary to investigate the 
coherence properties of condensates in terms of their mo- 
mentum distribution n(p), defined as the Fourier trans- 
form of the one-body density matrix pi{x,x') [3,S[3' 



n{p) 



dxdx exp ( — ip{x — x')) pi{x,x') , (22) 



where 



(23) 



Here the operators '^/{x) and ^{x.y - obeying the stan- 
dard bosonic commutation rules - annihilate and create, 
respectively, a boson at the point x; the average (...) is 
the ground-state average. Following Refs. it is 

possible to show that the momentum distribution n{p) 
can be written as 



n{p) — nQ{p)i 1 + a cos (pd) 



(24) 



Here uq^p) is the momentum distribution in the fully 
incoherent regime (which depends on the shape of the 
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FIG. 5; (Color online). Fisher information F of the ground- 
state \Eo) of the bosonic junction as a function of the pa- 
rameter C, = U / J for A'" = 30. Left panel: attractive bosons 
{U < 0). Right panel: repulsive bosons ([/ > 0). Dashed 
line: F obtained from the diagonalization of the Hamiltonian 
(HJ. Dot-dashed line: F obtained from Eq. HZ]) {U < 0) 
and Eq. (fT9)l (?7 > 0). The Fisher information F and C, are 
adimensional quantities. 




FIG. 6: (Color online). Coherence visibihty a of the ground- 
state |_Eo) of the bosonic junction as a function of the param- 
eter C, = U/J. Left panel: attractive bosons ([/ < 0). Right 
panel: repulsive bosons {U > 0). Solid hne: A'' = 20. Dashed 
line: A'^ = 30. Note that the horizontal axis of the right panel 
is in logarithmic scale. The coherence visibility a and are 
adimensional quantities. 



double-well potential Vdw{x)), and d is the distance be- 
tween the two minima of Vowix). a is a real quantity 
which measures the visibility of the interference fringes. 
This visibility is given by 



a 



N 



(25) 



and characterizes the degree of coherence between the 
two wells. 

The expectation value of the operator a^au is evalu- 
ated in the ground-state |£^o) and the visibility a is given 
by 



2 ^ 



,(o)„(o) 

■i H+ 



(26) 



We have analyzed a as a function of the parameter C, = 
U/J, for iV = 20 (sohd line) and TV = 30 (dashed line). 
The results of this analysis are reported in Fig. |6l As 
expected, we see that the coherent state, ^ = 0, has 
maximum visibility, a = 1. On the attractive side - left 
panel of Fig. |6]- it exists a value of (, say (a, for which 
the second derivative of a is zero. We have verified that 
Ca = C-Fi for a given TV such a quantity is close to the 
aforementioned (cr- We observe that both for iV = 20 
and for = 30, the coherence visibility a exhibits a 
plateau where it is almost independent on the interaction 
strength and very close to one. Except for this region, 
in correspondence of a fixed value of C, the greater is A^ 
the smaller is a. This behavior is reversed for repulsive 
interactions, as shown in the right panel of Fig. [6l from 
which we can sec that by increasing C, a slowly decreases. 
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FIG. 7: (Color online). Coherence visibility a of the ground- 
state \Eo) of the bosonic junction as a function of the pa- 
rameter ( — U/ J for A^ = 30. Left panel: attractive bosons 
{U < 0). Right panel: repulsive bosons {U > 0). Dashed 
line: a obtained from the diagonalization of the Hamiltonian 
Q. Dot-dashed line: a from Eq. ^ {U < 0) and Eq. (O 
{U > 0). The coherence visibility a and ^ are adimensional 
quantities. 

The Hellmann-Feynman theorem provides a relation 
between the coherence visibility a, and the derivative of 
the ground-state energy Eq with respect to J: 



leading to 



1 dEo 



N dJ 



(27) 



(28) 



By use of the perturbative results previously obtained for 
the ground-state energy iJo - see Eqs. P^ . ([T5)) . (PU)) - 
we can evaluate the coherence visibility ([28]) in the three 
regimes U/J^ — oo, U / J ^ +oo, and U ^ Q. In the 
first case, Eq. (|T6)) gives 



(29) 



(TV-i)C 



which, in the limit ^ —oo, tends to 0. For U/ J 
Eq. ([T8| provides 

A^-f 2 



-oo, 



(30) 
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FIG. 8: (Color online). Entanglement entropy S of the 
ground-state \Eo) of the bosonic junction as a function of the 
parameter ( = U/J. Left panel: attractive bosons {U < 0). 
Right panel: repulsive bosons ([/ > 0). Solid line: A'' = 20. 
Dashed line: — 30. Note that the horizontal axis of the 
right panel is in logarithmic scale. The entanglement entropy 
5 and C are adimensional quantities. 



which tends to zero when ( — > +oo. Finally, at weak 
coupling [U 0) the perturbative result (|20p gives the 
limiting value 



a 



1 



(31) 



In Fig. ([7]) we show the comparison between the coher- 
ence visibility a obtained by diagonalizing the Hamilto- 
nian ([T|). dashed line, and the ones obtained by using 
Eq. (|29p . for attractive bosons, or Eq. (PU)) . for repulsive 
bosons (see the dot-dashed line). 

At this point, we can conclude that the joint applica- 
tion of the Hellmann-Feynman theorem and of the per- 
turbations theory to calculate the ground-state energy 
makes possible to obtain very simple analytical formu- 
las for the Fisher Information F, Eqs. ^T7\ and p9)) . and 
for the coherence visibility a, Eqs. ([^^ and (|30p . at least 
over a certain range of values of the interaction strength. 

We analyze, now, the quantum entanglement pertain- 
ing to the ground-state \Eo)- When the system is in \Eo), 
the density matrix p is 



In this case, an excellent measure of the entanglement 
between the two wells is provided by the entanglement 
entropy S [sO]. This quantity is defined as the von Neu- 
mann entropy of the reduced density matrix defined 

by 



PL = Trjip 



(33) 



that is a matrix obtained by partial tracing the total 
density matrix ([5^ over the degrees of freedom of the 
right well. The entanglement entropy S is given by 



-TrpL log2 p'l 



N 

E 

i=0 



14"^ P log 



„(0)|2 



(34) 



For a given number of bosons N, the maximum value 
for S, say S(^jnax,N}, is given by log2(Af + 1). We have 
studied as a function of the parameter = U/J, for 
TV = 20 (solid hne) and N = 30 (dashed line). The 
results of this analysis - which is one of the most cru- 
cial studies of this paper - are reported in Fig. [8] We 
can see that the maximum entanglement is attained for 
a finite attractive interaction Cs- The larger is N the 
weaker is the attraction for which S has a maximum, 
which is equal to 4.26 when iV = 20 and to 4.72 when 



= 30. These values are smaller than S, 



[max ,20) 



4.39 



P^\Eo){Eo\ 



(32) 



and S(^rnax.30) = 4.95, respectively, and greater than 1, 
i.e. the entanglement entropy of a macroscopic cat state 
(sec Tab. 1). The macroscopic cat state is not the max- 
imally entangled state achievable in the bosonic Joseph- 
son junction. Moreover, we find that the interactions 
corresponding to the entanglement entropy maxima 
are (s - -0.13 for iV = 20 and Cs ^ -0.084 for N = 30. 
These values are close to the Ccr - onset of the coher- 
ence loss and the self-trapping regime within the junc- 
tion dynamics - discussed for the attractive ground-state 
coefficients (see the discussion about Fig. [5]). 

From the right panel of Fig. [5] we see that for repulsive 
interactions, S decreases by increasing (. 

We conclude this section by summarizing - for finite 
even number TV of bosons and at zero temperature - 
the expansion coefficients cf \ the Fisher information F, 
the entanglement entropy S, and the visibility a for the 
atomic coherent state ([3]), the separable Fock state ([1]), 
and the "macroscopic cat state" (O, see Tab. 1. 



State 



\ACS) 

\FOCK) 
\CAT) 



m 



2^ i\{N-i)\ 



'Ji,N/2 



F 



1 

TV 

1 



s 



loa 



->N 



1 ^ 



TV! 



TV! 2^-1 ^ i\{N-iy. 


1 



lOgn il 



Table 1. The coefficients 



the zero-temperature Fisher information F, entanglement entropy S, and visibility a for 
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the atomic coherent state ((3]), the separable Fock state Q, 
and the "macroscopic cat state" (O in the presence of a finite 
even number A'' of bosons. 



Results with two bosons 



When iV = 2, it is straightforward to calculate Ej and 



\Ej). The ground-state \Eq) reads 



\En) = 



)=a(^|2,0) + ^±^J+^|1,1) + |0,2)) , (35) 

where the A is the constant of the normalization which 
is 



A 



(36) 



The energy Eq is 



En = 



U - VU^ + 16 J2 



(37) 



The Fisher information p3|) calculated with respect to 
the state ([55]) is given by the following expression 



(38) 



16+ c + c^e + 16 



The coherence visibility ([25]) for the state ([55)1 is, instead, 
given by 



4(c + ^/eTiG) 



16 + + eve + 16 



(39) 



It is possible, also, to evaluate the entanglement entropy 
(p4|l when the system is in the state (j37| . We have that 
S can be written as follows 



S = -AH 21o 



2(C+VX^+16)' 



(16 + C2 + CVC^+16)' 



,C' + C\/C^+T6. 



l0g2 



(C + VC^TTe)^ 



32 + 2C2 + 2CVC2 + 16 



When the limit C ^ is met, the state (|35|) tends to 



V2|l,l) + |0,2) 



(42) 



which is the atomic coherent state ([3|) with two bosons. 
In this case, the energy ((37|) tends to —2 J. In the above 
mentioned limit, the Fisher information psp tends to 
1/2, the coherence visibility (|39p to 1, and the entangle- 
ment entropy (j40p to 3/2, in agreement with the first row 
of Tab. 1. 

Finally, let us analyze the limit C, +oo. We have 
that the ground-state ([SS]) tends to 



(43) 



l^o) = |l,l) , 



which is the separable Fock state (|4]) with two bosons. 
In the aforementioned limit, the Fisher information psp . 
the coherence visibility (P^)). and the entanglement en- 
tropy (go]) tend to 0. 



B. Quasi-classical coherent state with number 
fluctuations 

In this subsection, for the sake of completeness, we 
shall compare the TV-boson ground-state with the quasi- 
classical coherent state as follows. 

The familiar Josephson equations [l3| of the bosonic 
Josephson junction can be obtained from the Bose- 
Hubbard Hamiltonian ([IJ by using the quasi-classical 
coherent state \QC). This state is given by the tensor 
product of the coherent state \CS)l^ which describes the 
Bose-Einstein condensate in the left well, and the coher- 
ent state \CS) which describes the Bose-Einstein con- 
densate in the right well. The state ICS") is thus given 
by 



\QC) = \CS)l®\CS)r , 
where ICS")*.: (fc = L, R) is such that [U 



(44) 



(45) 



(40) 



As conclusive remarks, we observe that when ( —oo 
the state (l35ll tends to 



The complex quantity Zk is the eigenvalue of the annihi- 
lator in the kth well, i.e. 



dk\cs)k^z,\cs)k 



(46) 



V2 



2,0) + |0,2) 



(41) 



The absolute values of the two are related to the av- 
erage occupancy of the two wells: 



This is the macroscopic superposition state ([S]) with two 
bosons. In the limit C — oo, the energy tends 
4J2 

to [/ H . The Fisher information (1551) tends to 1, 

U 

the coherence visibility (P^)) to 0, and the entanglement 
entropy (gHl) to 1. 



Nk = {QC\hk\QC) = \zk\ 



(47) 



So 2fc are conveniently parametrized as Zk = 
\/Nf;^ exp(i0fc), where Ok are phase variables. It is im- 
portant to stress that the quasi-classical coherent state 
\QC) does not have a fixed number of particles and 
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0.05-0.025 



FIG. 9: (Color online). On the vertical axis, the energy of the 
ground-state |-Bo), -Bo (solid line), the energy (|48|) calculated 
for z — 0, E.^o (dashed line), the energy (|48p calculated for 
z = Zsb, -Ez=Zj,5 (dot-dashed line). On the horizontal axis the 
parameter ^ — U/J. Left part of the upper panel: attractive 
bosons {U < 0). Right part of the upper panel: repulsive 
bosons (U > 0). Lower panel: three insets for attractive 
bosons. We have set A'' = 30. E is in units of J and is an 
adimensional quantity. 



consequently it allows for non vanishing quantum fluc- 
tuations of the total number of bosons. The average 
boson number in the state \QC) is simply given by 

7V = iVi+iVfl = |zi|2 4-|zfl|2. 

The expectation value E of the BH Hamiltonian ([1]) 
evaluated with respect to the state (03]) is 



E = {QC\H\QC) = -NJ^/T 



(1 + -^) 



(48) 

where z = {Nl - Nr)/N and (f> = Or - Ol- In the 
formula (|48p we recognize the total energy of a nonrigid 
pendulum, which represents the mechanic analog of a 
bosonic Josephson junction d, [l^. In fact, the state 
\QC) is expected to be close to the ground state \Eq) in 
the weak-coupling regime jl7| . It is straightforward to 
show that the minima of the energy are obtained 



with (/) = and 

Z = Zsb = 

z = if C> 



(2 JSf2 

2 

'n 



if c< 



N 



(49) 

Note that C = —2/N signals the onset of the self- 
trapping regime within an attractive ID bosonic Joseph- 
son junction [l^. The energy pS)) calculated for z ~ 0, 
say i?z=Oj 3'iid for z = z^i,, say Ez=zsb^ becomes 



E, 



E, 



'NU 
N\ — 1 



U 



(50) 



(51) 



It is interesting to compare the energy E of the quasi- 
classical state \QC) with the energy Eq of the state |i?o), 
that is the ground-state of the BH Hamiltonian with a 
fixed number N of bosons (see also Rcf. [13 )• ^From 
the upper right panel of Fig. [9] - obtained with N — 
30 - we see that for repulsive bosons, £^o (solid line) is 
smaller than Ez=q (dashed line). For attractive bosons, 
the situation is much more complex. Let us focus on the 
lower panels of Fig. [9l We can see that there are three 
different regions. In the first one (left panel: C, < ~2/N~ 
A (A > 0) ) we have that Ez=z,t < Eq < Ez=o; in the 
second region (middle panel: ~2/N ~ A < ( < ~2/N ) 
we have that Ez^z^t < Ez=o < Eq; in the third region 
(right panel: C > ~2/A^) we have that Ez=o < Eq. Note 
that the greater is N the smaller is A. 
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FIG. 10: (Color online). Coherence visibility a of the bosonic 
junction as a function of the dimensionless parameter ^ = 
U/J. Left panel: attractive bosons {U < 0). Right panel: 
repulsive bosons {U > 0). Solid line: results from Eq. (|26[) . 
Dashed line: results from Eq. (|52|l . Note that the horizontal 
axis of the right panel is in logarithmic scale. The coherence 
visibility a and ( are adimensional quantities. 

The coherence visibility calculated with respect to 
the state \QC) reads 



a 



(52) 
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with z given by Eqs. (|49p . In Fig. [10] we report the 
comparisons between the coherence visibiUty calculated 
by using the ground-state |i?o) of the BH Hamiltonian 
([1]), and that one calculated by using Eq. ([5^. Eq. 
([52)1 . Fig. [To] clearly shows that the coherence visibility 
a of the quasi-classical coherent state \QC) is strictly 
equal to one only for a non-negative interaction strength. 
Moreover, from the left panels (attractive bosons) of Fig. 
[TOl we can see that the greater is the number of particles 
in the system the better is the agreement between the 
predictions deriving from the two approaches. 

The Fisher information F, Eqs. p^ - (|T^ . calculated 
with respect to the state |QC), is simply given by 



This is the same result that one obtains when the sys- 
tem is in the atomic coherent state ((31), i.e. the exact 
eigenstatc for [/ = 0. 

V. CONCLUSIONS 

We have investigated an atomic Josephson junction 
in the presence of a finite number of interacting ultra- 
cold bosons confined in a quasi one-dimensional geom- 
etry. By using the two-site Bose-Hubbard Hamiltonian 



toolbox, we have carried out the zero-temperature anal- 
ysis by finding the ground-state of the system. We have 
characterized the presence of macroscopic Schrodinger- 
cat states, the coherence properties, and the quantum 
correlations by calculating the Fisher information, the 
visibility of the interference fringes in the momentum dis- 
tribution, and the entanglement entropy as functions of 
the boson-boson interaction strength. The joint applica- 
tion of the Hellmann-Feynman theorem and of the per- 
turbative approach to evaluate the ground-state energy 
allows to obtain analytical formulas - within the interac- 
tions range in which the perturbative theory is applicable 
- both for the Fisher information and for the coherence 
visibilty. The Fisher information and the coherence visi- 
bility predicted by these formulas are in good agreement 
with those evaluated by diagonalizing the two-site Bose- 
Hubbard Hamiltonian. In the attractive regime we have 
found that the presence of macroscopic Schrodinger-cat 
states corresponds to a maximum of the Fisher infor- 
mation, while the entanglement entropy has a maximum 
when the system begins to lose coherence and to exhibit 
self-trapped dynamics. The maximal entanglement en- 
tropy is larger than that of a macroscopic cat state and 
in fact very close to its theoretical upper limit. 

GM thanks Alessio Serafini and S.M. Giampaolo for 
useful discussions. 
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